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We study the normal state and the superconducting transition in the Attractive Hubbard Model in 
three dimensions, using self-consistent diagrammatics. Our results for the self-consistent T-matrix 
approximation are consistent with 3D-XY power-law critical scaling and finite-size scaling. This is 
in contrast to the exponential 2D-XY scaling the method was able to capture in our previous 2D 
calculation. We find the 3D transition temperature at quarter-filling and U = —At to be Tc = 0.207i. 
The 3D critical regime is much narrower than in 2D and the ratio of the mean-field transition to Tc is 
about 5 times smaller than in 2D. We also find that, for the parameters we consider, the pseudogap 
regime in 3D (as in 2D) coincides with the critical scaling regime. 
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Inspired by the discovery of high-temperature super- 
conductors there has been considerable interest for some 
time in understanding how the BCS picture gets modified 
in the case of short coherence length superconductivity. 
A simple model for a correlated superconductor is the 
Attractive Hubbard Model that allows one to tune the 
attraction to an intermediate regime between the BCS 
and Bose limits where the fermions are still degen- 
erate but the coherence length starts to approach the 
interparticle spacing. Earlier studies ^ of this model in 
2D using Monte Carlo methods revealed the development 
of a pseudogap in the normal state. In this paper we em- 
ploy self-consistent diagrammatics that is a variant of 
the Fluctuation-Exchange Approximation ||-||]. Specif- 
ically, we use the self-consistent T-matrix approxima- 
tion (STA) that only retains fluctuations in the particle- 
particle channel. These are presumably the only relevant 
fluctuations in the model for the densities we consider. 
We have found before that quite remarkably, the STA 
approach in 2D captures the critical 2D-XY fluctuations 
with a Kosterlitz-Thouless (KT) transition temperature 
in good agreement with Monte Carlo Q . In addition we 
found a very broad critical regime that coincides with 
the development of a normal state pseudogap. It is very 
surprising that self-consistent diagrammatics can capture 
critical behaviour relatively well. In order to understand 
this better, and encouraged by the success in 2D, we 
now study this model in three dimensions where tech- 
niques like Monte Carlo are severely limited by system 
size constraints. 

We find that the STA method indeed yields results in 
three dimensions that are consistent with the expected 
power-law scaling rather than the exponential KT-scaling 
we found in 2D. In addition, the success of finite-size scal- 
ing allows us to determine Tc in 3D. In addition to observ- 
ing a much narrower pseudogap regime than in 2D, we 
can quantify the difference in pairbreaking fluctuations 
by looking at the ratio of the mean-field (1-loop) Tc to 
the actual transition temperature and find that this ratio 
drops from 12 to about 2.6 in going from 2D to 3D, for 
the parameters we consider. We also conclude that since 
the dramatic KT-scaling observed in 2D is not seen in 



3D, this result is not a numerical artifact. 

While superconductivity in the Attractive Hubbard 
Model is of course s-wave, it is fair to assume that the 
strong correlation effects in this model will also apply to 
a more realistic c?-wave theory. In fact, an advantage of 
the STA is that one can explicitly consider models with 
(i-wave attraction In the context of the cuprates 
our work suggests that the loss of low-energy spectral 
weight in the pseudogap normal state can be attributed 
to pre-cursor pairing as observed in our calculation. By 
studying dimensionality, the observation then is that un- 
derdoped cuprates behave more like our 2D calculations 
while optimally and overdoped materials have pseudogap 
behaviour more in line with our 3D results. Since, within 
a given family of materials, the anisotropy can increase 
several hundred-fold upon underdoping, the implication 
is that this change in effective dimensionality plays an 
important role in enhancing the underdoped pseudogap. 

Our starting point is the hamiltonian for the Attractive 
Hubbard Model 



hc-\U\J2', 



(1) 



where the tight-binding hopping sums over nearest neigh- 
bour pairs on a three-dimensional square lattice and the 
potential captures singlet pairing for a short-range po- 
tential. The standard perturbative approach is to eval- 
uate the pair susceptibility at the 1-loop level. Such a 
Gaussian theory amounts to re-summing the BCS ladder 
diagrams yielding an RPA-likc form 



Ho(q, iujri) 

l-\U\UQ{q,iuJn) 



(2) 



where iujn are Matsubara frequencies and Hq is the bare 
bubble in terms of free fermion propagators. In the 
limit of small U , this framework reduces to BCS-theory. 
We are interested in intermediate attraction where the 
fermions are still degenerate but the pairing is strong. 
The value U — —At that we use throughout this paper 
is representative of such intermediate attraction. The 
model also deviates strongly from BCS theory near half- 
filling, where the Attractive Hubbard Model maps into 
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the Hubbard Model with on-site repulsion. Since we do 
not wish to include the effects of the associated CDW 
instability at half-filling, all our results in this paper will 
focus on quarter-filling. 

The superconducting instability signaled by the Thou- 
less criterion that determines the transition tem- 
perature as the highest temperature that satisfies 
l/XgQus(r = Tc) = 0. 

Our approach is to extend this framework to the self- 
consistent T-matrix (STA) formulation, which simply 
generalises the bare bubble in (j^) with dressed propa- 
gators for the intermediate fcrmion states: 



n(q, iun) 



^ G{q-k,iuJn-ikn)G{k,ikn) (3) 



where the dressed Green's function G(k, ifc„) = [ikn — 
Ek + — ikn)]"^ is defined in terms of the pair sus- 
ceptibility itself, through the self-energy: 

E(k,zfc„) = 4 X! G{c[-k,iuJn-ik„)T{q,iujn)- (4) 



f3 

Technically, we explicitly keep only corrections of order 
and higher which implies that the Hartree term is 
already included in the unperturbed chemical potential. 
This higher order contribution is expressed though the 
effective interaction: 



r(q, iun) 



n(q, lUJn) 
1- |i7|n(q,za;„)' 



(5) 



The self-consistent solution of equations (H) and (Q) 
are achieved numerically, by considering a finite A^x A^x N 
lattice with discrete momenta and only retaining a fi- 
nite number of Matsubara frequencies jlOj . We start off 
with the unperturbed solution and iteratively solve these 
equations (which are convolutions) using FFT's until self- 
consistent convergence is established. 

Perturbatively this of course amounts to re-summing 
an infinite series of diagrams. At first one may expect 
that the self-consistent feedback will broaden the inter- 
mediate fermion states in and simply re-scale values 
for physical quantities resulting from a 1-loop calculation. 
This is in fact not the case as can be seen from the results 
in Fig. 1 for the inverse of the uniform static pair sus- 
ceptibility x('Z = 0, iti'n = 0) evaluated for quarter-filling 
{n-ia) = J and the intermediate attraction U = —At. The 
results are for a 32x32x32 lattice. On the RHS the 1- 
loop resuh intersects the axis at TJ"'"'*(3i:') = 0.532t. 
This result is essentially independent of lattice size and 
the 1-loop approximation yields a non-zero transition 
temperature even for finite systems. In any exact treat- 
ment, one of course only expects phase transitions in the 
thermodynamic limit. We also indicate the 2D transition 
temperature in this 1-loop approximation for the same 
parameters at TJ"'"^(2L>) = 0.6111. 



0.20 
0.15 
g0.10 
0.05 
0.00 
-0.05 







( ) 





T ' ' ' — n ' ' 

-*;:3DSTA ,' 


- 




) 




-■ 3D gaus i 


- 


( 


) 




; 




( 

( ) 










i 






; — ^ 


TJ2D) 








;' Tr(2D) 



0.0 



0.2 



0.4 
T 



0.6 



FIG. 1. Inverse pair susceptibility for a 3D lattice of linear 
length L = 32, density (nia) = j and U = —it. The symbols 
indicate the STA results and the broken line on the right is 
for the 1-loop (mean-field) theory. Also shown are the 2D 
mean-field and KT transition temperatures. 

The 3D STA results indicated by the circles show a 
fundamentally different behaviour. The results for the 
inverse susceptibility never cross zero - indicating that 
the self-consistent calculation reflects the fact that there 
is no transition for a finite system. This is of course 
what one would expect from an exact treatment - but 
is missed by 1-loop theory. Simply put, the self-energy 
corrections broaden the intermediate fermion propaga- 
tors sufficiently to avoid the superconducting instability. 
Thus the STA is fundamentally supcriour to 1-loop in 
that it captures the physical requirement that there are 
no phase transition in finite systems. 
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FIG. 2. Finite-size dependence of x{T, L) in 3D for 
{■n-i^) = J and U = -U. 

Is there any more to this result? Further investigation 
reveals that the low-temperature static uniform pair sus- 
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ceptibility xC^i ^) for the STA depends on system size as 
shown in Fig. 2. Just as one would expect from funda- 
mental principles, the susceptibility increases towards a 
maximum that is controlled by system size and can only 
diverge in the thermodynamic limit. Since the formalism 
does not include anomalous propagators, it is limited to 
the disordered side of the phase diagram. This plot is 
quite similar to our earlier results in 2D, except that 
the attempt at a transition appears sharper and is at a 
higher temperature in 3D. 

From this finite-size data we use standard scaling ar- 
guments | [Tl|] to extract the transition temperature in 
3D. For a finite system there are two length scales 
namely the correlation length ^ and the system size 
L. One expects from a 3D scaling hypothesis that 
X(T, L) = \T-T,r fi{L/aT)) where ^{T) is the cor- 
relation length in an infinite system. From ^(T) ~ 
\T — Tc\~'^ and the scaling law 7 = i/(2 — 77), it follows 
that the following product is only a function of the ratio 
of the two length scales: 



x(r,L) = /(i/e(T)). 



(6) 



At Tc the correlation length diverges and the argu- 
ment on the RHS becomes zero, independent of system 
size. Plots for this function for our STA data on different 
system sizes are given in Fig. 3. The numerical data is 
consistent with the finite-size scaling hypothesis and we 
can estimate Tc ~ 0.207i for our standard set of param- 
eters. This value of T is indicated by the vertical line 
where the plots become independent of L. 
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FIG. 3. Finite-size scaling of f{L/^{T)) in 3D for 
(jiio-) = J and U — —4t. The common intersection yields 
n = 0.207t. 



scaling regime that coincides with the regime in which a 
pseudogap develops in the density of states at the Fermi 
level, No{T). We also found very good finite-size scaling 
in our data. Given the difficulty of theoretically captur- 
ing the correct critical behaviour this was very surprising 
but one could always question whether this result was a 
numerical artifact. 

Here we have demonstrated, that in 3D the STA data 
is again consistent with a finite-size scaling hypothesis - 
yielding the first estimate of Tc = 0.207t for the param- 
eters in question, in a calculation that goes beyond the 
mean-field or 1-loop class of approaches. 

It is interesting to note that while the 1-loop estimate 
for the transition temperature is higher in 2D than in 
3D for the same parameters (due to density of states 
effects) the opposite holds for the actual transition tem- 
peratures. We conclude that the enhanced fluctuations 
in 2D overcome the density of states effect. If we con- 
sider the ratio of mean-field to actual transition temper- 
atures R = Ta^-^'/Tc we find R2D = 0.611/0.05 - 12 vs. 
RsD = 0.532/0.207 ~ 2.6. 

Let us next consider the bulk scaling behaviour. Is this 
consistent with the expected 3D-XY result 
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\T-Tc 



(7) 



with 7 ~ 1.33? In 3D we are limited to smaller sys- 
tems than in 2D and the transition itself is narrower. 
This limits the accuracy with which we can determine 
the critical scaling in the thermodynamic limit within 
our calculations. In Fig. 4 we plot the finite-size data 
X{T, L) against the RHS of (g). 
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FIG. 4. Finite-size results for x{T, L) and the expected 
scaling form in (mj with exponent 7 ~ 1.33. 



In our earlier 2D calculation, the STA yielded results 
that were very consistent with the Kosterlitz-Thouless 
scaling of the 2D-XY model, with an estimate of T^t — 
0.05t in good agreement with previous Monte Carlo es- 
timates [pi. In addition, we found very broad critical 



For each system size, the high-temperature data coincides 
with the scaling form in and at lower temperature one 
observes the expected finite-size corrections. While the 
data is clearly consistent with power-law scaling and the 
fit with the exponent 7 ~ 1.33 does better than the mean- 
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field 7m _F = 1, our accuracy is insufficient to determine 
the exponent a-priori. In addition we can estimate the 
width of the critical regime as extending about O.lt above 
Tc, where the high-temperature results start to deviate 
from the asymptotic form. 

While we cannot determine the power-law exponent, 
the 3D data is dramatically different from our previous 
2D result that found data consistent with the exponential 
KT scaling expected in 2D. This gives some confidence 
that the 2D result was real and that the STA method 
captures the expected critical behavour in both 2D and 
3D quite well. This implies that for this particular model, 
self-consistency introduces scaling corrections that go be- 
yond dimensional analysis - a feat that normally requires 
a technique such as the renormalization group. We be- 
lieve that this may be connected to the rather simple 
C/(l) symmetry that superconductivity breaks, and that 
this result does not apply to more complicated problems. 

Let us finally turn to the normal state pseudogap. The 
density of single-electron states can be obtained from the 
Green's functions dressed with our calculated self-energy 
corrections S(k, ikn) and continued to retarded frequen- 
cies using Pade analytic continuation. In Fig. 5 we show 
the temperature-dependence of the density of states at 
the Fermi level for our standard density and interaction 
strength. We show results for both 3D and 2D and in- 
dicate the respective T^s by the vertical lines. One can 
clearly see that the 3D pseudogap is significantly smaller 
than in 2D and again coincides with the critical regime. 
In earlier analytic work we calculated the 3D pseu- 
dogap for a continuum model within the 1-loop theory 
and found that even at this level, the pseudogap correc- 
tions go beyond the expected results in a small energy, 
long-wavelength expansion. 
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haviour for the 3D-X'K model. At quarter-filling and 
for and intermediate attraction U = — 4t, we find in 3D 
Tc ~ 0.207i, about 2.6 times smaller than the mean-field 
transition temperature, and a narrow critical and pseu- 
dogap regime. We also find that while the mean-field 
transition temperature is higher in 2D than in 3D, the 
opposite holds for the actual Tc. Our results suggest 
that the previous 2D results are not a numerical artifact 
and that the STA formalism can distinguish between the 
fundamentally different nature of fiuctuations in two and 
three dimensions. 

We wish to acknowledge useful discussions with 
Alexander Nazarenko who contributed in a key way to 
our earlier work using the STA technique. 
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FIG. 5. Pseudogap in the density of states at at the Fermi 
level in 3D and in 2D for (riia-) = j and U = —4t. 



In summary, we have shown that results for the self- 
consistent T-matrix approximation in 3D are consistent 
with the expected finite-size and power-law critical be- 
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